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Abstract. We present a new theoretical approach to describe the quantum behavior 
of a macroscopic system interacting with an external irradiation field, close to the 
resonant condition. Here we consider the extremely undcrdamped regime for a system 
described by a double well potential. The theory includes both: transitions from one 
well to the other and relaxation processes. We simulate resonant phenomena in a 
rf-SQUID, whose parameters lie in the range typically used in the experiments. The 
dependence of the transition probability W on the external drive of the system ipx shows 
three resonance peaks. One peak is connected with the resonant tunneling and the two 
others with the resonant pumping. The relative position of the two peaks correlated 
to the resonant pumping depends on the pumping frequency v and on the system 
parameters. The preliminary measurements on our devices show a low dissipation 
level and assure that they are good candidates in order to realize new experiments on 
the resonant phenomena in the presence of an external microwave irradiation of proper 
frequency. 
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1. Introduction 

Josephson devices are convenient instruments for investigating macroscopic quantum 
effects [1]. First experiments focused on incoherent plienomena such as macroscopic 
quantum tunnehng (MQT) [21 [U HI [5] and energy level quantization (ELQ) [SI [3, [8]. 
Recent experiments have shown a coherent superposition of distinct macroscopic 
quantum states in SQUID systems [9l[T0l[TT] under the action of an external microwave 
irradiation. These results have stimulated many researchers to search for new quantum 
phenomena at macroscopic level as well as a method to distinguish them from classical 
predictions. Moreover quantum effects in SQUID dynamics are interesting because of 
possible apphcations to quantum computing [12] , as demonstrated by recent experiments 
[To] , [11] . Anyway these results need a final confirmation, since some of them have been 
reproduced by using a classical model to describe the phenomena |13j . 
The main problem to observe quantum phenomena at macroscopic level is the interaction 
with the environment which produces dissipation of coherent effects into classical 
statistical mixture of quantum states. Among a great variety of incoherent effects, 
the environment induces a finite width 7 of quantum levels [H] and, in this situation, 
resonant phenomena can occur only if 7 is small compared to the energy difference 
between levels, h'j <C hflp, where 



is the plasma frequency, and v?mm is the coordinate corresponding to the minimum of 
the considered potential well. We have already discussed [T5| the behavior of microwave 
induced transitions in a rf-SQUID in the moderate underdamped regime, when 7 is larger 
than the probability of quantum tunneling under the potential barrier, T/v, (7 ^ T/v) 



In the present paper we extend the theoretical analysis to the small viscosity limit, 
corresponding to 7 ~ T/v, referred as the extremely underdamped regime. In this 
conditions we can neglect the interwell transitions with respect to the tunneling 
transitions. Up to now, this regime has been investigated only for other Josephson 
devices [TOl [8]. Here we present a new theory specific for the rf-SQUID, characterized 
by a double well potential (FigH]), in order to identify coherent phenomena and to 
suggest a way of measuring them. 

It's worth nothing that Averin et al. [17], [18], [19] already reported a theoretical 
investigation of resonant quantum phenomena in rf SQUIDs. Nevertheless in that 
paper all quantities are postulated and not found from the initial Hamiltonian, while 
we don't not use any fitting parameters, but only external parameter values entering in 
the Hamiltonian. 

We study the interaction between a macroscopic quantum system (the rf-SQUID) and 

an external microwave irradiation, for frequencies close to resonant conditions. 

The probability of a transition under the potential barrier decreases exponentially with 




(1) 
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Figure 1. A sketch of an rf-SQUID double-well potential U((/9), with two minima 
located at 1^9™*" and (^s™™- 'Ptop is the position of the potential maximum. The 
coordinates 1^1,1^2,^3,'^ a are the "turning points" of energy Ef-^. In the same way 
it's possible to define turning points for energies Ef^, Eq, Er, El. We study the case 
E]^ = Eq, that is the ground state in the left potential well. 

the decreasing energy of the quantum state. Tunnehng occurs from a level close to the 
barrier top. In this case, as we'll show in the following, the transition probability curve 
can present three peaks. 

A further requirement to observe the quantum tunneling phenomenon is that the number 
of levels n is not too large, so that no crowding of levels happens. At the same time we 
assume the number of levels in the potential well is large enough, so that the motion 
of a particle with energy close to the top of the potential barrier can be considered 
quasiclassical. In this limit we use approximations reported in |20j . 
In the quasiclassical limit, the density of levels close to the barrier top increases, but in 
the considered case this effect is not significant, since numerical factors are of the order 
of 2^ Inn [20]. Furthermore, near to the barrier top, repulsion of levels is relevant and 
hence the equation for the system wave functions should be solved exactly. Finally the 
tunneling probability of a particle through the potential barrier depends very strongly 
on the particle energy E and on the energy difference between neighboring levels. The 
position and the width of each level depends on the external parameters of the system. 
A sketch of the double well potential of the rf-SQUID here studied is reported in FiglH 

Each macroscopic flux state of the rf-SQUID corresponds to the system wave function 
confined in one of the two potential wells. By varying the external parameter ip^ it's 
possible to control the energy difference between the levels in different wells. 
Out of resonance conditions, the energy levels are localized in each well. When two levels 
in different wells are in resonance, a coherent superposition of distinct states occurs and 
a gap appears between the energy levels. As a consequence of the coherent quantum 
tunneling process, in exact resonant conditions, the two levels have an energy tunnel 
splitting A, and the wave functions spread over the two wells, generating delocalized 
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Figure 2. The four essential energy levels considered in this study. All the quantities 
are referred to the barrier top. Here Ef-^ and Ef^ are the energies of the delocalized 
levels. £^0 is the first level in the left potential well below Ef^ and En is the first level 
in the right potential well below Ef^. 



states with energies Ef-^ and Ef^, as shown in Fig|2l It's worth noting that the 
energy tunnel sphtting A corresponds to the tunnehng probabihty T/v in exact resonant 
conditions and can be experimentally observed only if resonant levels are very close to 
the barrier top and in small viscosity limit |9j. 

Moreover, to study the tunneling process from the upper levels, it's necessary to populate 
them through a pumping process from the ground state of the left well: this requires 
the application of a microwave irradiation of proper frequency u = uj/2tt and amplitude 
I. Even a small amplitude of external microwave irradiation will cause an essential 
change in the population of excited states and will therefore produce a large effect in 
the tunneling probability. 

We study the rf-SQUID dynamics in these particular conditions and we show some 
numerical simulations. For a fixed value of the pumping frequency u, the transition 
probability curve W from one potential well to the other one (as a function of the 
external parameter ip^) can present three maxima: one maximum is connected with 
the resonant tunneling and the two others with the resonant pumping toward the two 
non-localized levels with energies Ef-^, Ef^. 

In the considered phenomenon, friction, that is due to the interaction with the external 
environment, has a relevant role since it leads to a finite width 7 of levels and destroys 
the coherence [21]. In the following sections we will discuss the induced resonance 
phenomena by using the density matrix formalism [22] . 

In particular in section II we discuss the transition probability function in presence of 
resonant pumping, in section III we show how to calculate levels position and transition 
matrix elements, in section IV we determine the energy spectrum in the vicinity of the 
crossing point and , finally, in section V we show the results of numerical simulations 
and first experimental data for the dissipation level for these rf-SQUID based devices. 
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2. Transition probability function in presence of resonant pumping 



In the zero approximation of the system interaction with the thermal bath, the 
hamiltonian describing the rf-SQUID in the presence of an external microwave with 
amplitude X is 



n 1 / N 

+ Pl cos cos{ujt)(p 



(2) 



2M dif^ 

where X has the dimensionality of a current, M is the "mass" of the Josephson junction 
defined as 



M 



(3) 



and C is the junction capacitance. Following a common practice|l], the dimensionless 
variable (p^ is obtained by referring the magnetic fiux to $o/2 and normalizing it to 
$o/2vr. In eq.(I2]) quantities Uq and Pl are free parameters of the system and are defined 



as 



Uo 



$0 

277 

27rLL 



1 
L 







(4) 



(5) 



L and Ic are the inductance and the critical current of the rf-SQUID, respectively. The 
simplest case is when Uq and Pl are constant, while ipx changes. 

Our considerations are valid for any potential U{(f). The essential point is that the 
potential U{ip) has the form given in FigH] in the proximity of the points ip. All the 
parameters (p^, Uq, Pl depends on the structure of the considered SQUID. 
Here we choose parameters such that, near to the top of the barrier, there are two 
close levels Ef^, Ej^, that can be considered as belonging to both potential wells 
simultaneously (delocalized levels) as shown in FiglH Moreover it's worth noting that 
the quantity is proportional to 

The dissipation is accounted for by an effective resistance -Re// in the RSJ model for the 
junction [23.]. So a bigger -Re// implies a better insulation, that means a smaller width 
of the levels. Here we consider an rf-SQUID whose Josephson junctions presents a large 
value of resistance -Re//, so that the condition 



h 



> 7i,72 



(6) 



is satisfied, where 71,72 are the widths of levels Ef^^, Ef^ respectively. 

To describe the resonant tunneling induced by the resonant pumping quantitatively, we 

use the following system of equations [H] for the density matrix elements pj- 



— = -cosM)^ 



-{m\ip\f)exp 



(j|v3|m)exp 



Em — Ei 



Em. — Ej 



mm 
n In 



Pf 
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Here j, m, f, n are the considered levels, Ej, Em, Ef, En are their energies and t is the 
time. 

The matrix elements Wj^ entering in eq.(I7]) are defined in Appendix A and C. They 
are nonvanishing provided that 

I {Em - Ej) - {En -Ef)\<^ hQp (8) 

If these conditions are satisfied and if the temperature is low enough, the nonvanishing 
terms Wj^!^ of the transition probability matrix are given by the following equation: 

W}: = ^^-^ (^1 +tanh|,^ Ar(cD)- [(j|e^|m)(/|e-^|n) + (jle"^ |m) (/|e^ (9) 
where 

^ ^ {Em - Ej +En- Ef) ^^^^ 

and N {uj) is defined as 

N {Co) = ^coth(^-^^ (11) 

Only the four non-diagonal elements of the density matrix pj^, p^^, Pq\ Pq^ are nonzero, 
and they are connected by the expressions 

P'o = {pIY p'o = {P%y (12) 

In presence of an external pumping also the two diagonal elements pj-^ and pj-^ are non- 
zero. 

We suppose that the temperature T is much smaller than the plasma frequency 
{ksT -C hflp), such that it is of the order of the energy difference ksT ^ {Ef-^ — Ef^). 
In the present article we consider T = 50 mK. When these hypotheses are satisfied, 
from eq.© we obtain 



dt 2e 
dpi 



cosM)(0|^|/i)er' VpO + wl%p%-^,pl (13) 



-- cos(^t)(0|^|/2)eV- ^ ;pO + Wl'f^pl - 72P;, (14) 



dt 2e 

where the widths 71, 72 of levels are given by the following expressions 



1 {Wj^ll + Will + Will) (15) 



^1 

The indexes L and R refer to the first energy levels below Ej^ in the left and right wells, 
respectively. The tunneling transition probability between two states in different wells 
decreases very quickly for states below the barrier top, such as e~^'^", where n is the 
number of states counted from the barrier top [23]. As a consequence we can neglect 
transitions into lower levels in the right potential well. 

In order to study this phenomenon, only five energy levels (-E/i, Ef^, Eq, E^, Er) are 
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essential quantities. Eq is the energy of the ground state in the left potential well, 
is the first level in the left potential well below Ef^ (in the simulations presented here 
El=Eo), Er is the first level in the right potential well below Ef^. Energy levels as a 
function of the external parameter (/jj,, are shown in Figj2l The pumping frequency v is 
supposed to be close to the frequencies Ef^/h and Ej^/h (referring both to Eq/K). In 
our simulations E^ is also the ground state of the left potential well (it's the same as 
-Eq), while in the right well there are many more levels, not involved in the described 
process. 

If these conditions are satisfied, the solutions of eqs. (fT3|) . (13) are: 



Ph = ^/i^ 
pI = Bf.e 



l\ UJ i\ \t 



Bf.e 



(16) 



+ Af,e 

In eqs. (fT6ll coefficients A/^, Af^, Bf^^, Bf^ are defined as follows 



4e 



%2 - 



h 



■ «72 - 



—n 



(17) 



5/2 



h 



- ni 



- «72 



The diagonal elements and p^^ of the density matrix can be find by solving the 
following equations 



— ^ = COs(co't) 

dt 2e ^ ' 



dt 2e 



— cos(ti;t) 



pWvW 

P/.(/2|^|0)e 



Po^MO|^|/i)e 



.Ef^-Eo 



+ WflflpH - Hl4^ 



These rate equations completely describe the dynamics of the process. The solutions of 
this system are respectively: 



pj; (t) = pg + T,e 

Pfl (^) = P/2 + ^le' 



/l 



(20) 



+ P*e V " ; (21) 

The expressions for J^i and Vi are given in Appendix A. By using equations (!20|) and 
(j2Tl) inside equations ( fTSl) and (fT9l) . we obtain the expressions of the quantities p^^ and 
pj^^ as following: 

t2 1 
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Finally we can introduce the transition probability W from the left to the right potential 
well, is a function of the external parameter cpx and computing this quantity is the 
main goal of the paper. The transition probability is given by the expression 

^ = + <iV/2 (24) 

where, in the low temperature limit, here considered, from eq.Q we know that 
Rfi _ 2 (% - Eji) 



W 



Rfi 



w 



Rf2 



Rf2 



ReffC^ 
ReffC^ 



{R I e^^/^ 
{R I e'^/^ 



/i) P 

/2) r 



(25) 
(26) 



while the expressions for p^^ and given in eqs. fl20l) and (pT]) respectively. In the 

eqs. ([22D and (12^ . the quantities JF^ and Pi determinate the value of the oscillating 
part of the tunneling probability W. 

It's worth noting that change of the external parameter (fx leads to a redistribution of 
wave functions relative to energies i?/^, Ef^ between left and right potential wells, so 
that transition matrix elements and widths of levels 71, 72 are functions of the external 
parameter (fx- To complete our considerations we will find the levels position and the 
transition matrix elements as a function of the external parameters, as shown in the 
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next section. 



3. Levels position and transition matrix elements. 

As we have seen before, a necessary condition for resolving the splitting is that the 
experimental linewidth of states (71, 72) be smaller than A (A ^ hji, ^172): it can be 
satisfied only for levels near to the barrier top (Fig{T]). In Fig{T] y^i, (y92, V^s, V^4 are the 
"turning points" , that are the solutions of equation 

t/(<^l,2,3,4) = E (27) 

where E is the energy value. In Figil] the points yj^*", yj^*" are coordinates 
corresponding to the minima of potential U{(p) and the point (ptop is the coordinate 
corresponding to the maximum of the potential U{ip). For energy values E close to the 
top of the barrier, the wave function can be written in the form 

1 




^ = ■ -77 sin - + / V2M (E - U ((p))d(p in the left potential well 

{2M{E-U{^))f^ ' ' ' ' 

^ = A- + i){2MUiY^^{ip - iptop)) + A*V_i^ ((1 - i){2MUi)^^^{ip - cptop)) near the barriei(g8]) 

B I TT 



^ = -77 sin - + / V2M (E - U((p))d(p in the right potential well 

In eqs.f l28l) Vp is the parabolic cylinder function [25], while quantities A and B are 
numerical factors, G is a normalization factor, defined in the following. The other 
quantities are defined as 

Utop = U{^top) (29) 

(30) 



2\dip 



-Utrm — E 



X = ^2MUi^^ (31) 



First expression in eqs. fl28|) is valid in the left potential well, the second one is valid 
close to the barrier top, the third one is valid in the right potential well. Coefficients 
A, B are defined by imposing boundary conditions for these three expressions. Near to 
the top of the potential barrier we use the exact solutions of the Schroedinger equation 
for the considered system, whereas in the vicinity of the points yji, (p^ we use the well 
note method of outgoing in complex plane [26]. The relation between coefficients A and 
B can be found by matching coefficients in the right potential well 



B 

23/4(2Mf/i)V8 ^""P 



ttA ivr i\ i\ ^ f 2\ . f , , ^ , ^ , 

— -- + - + -In (^-j dip^2M{E-U{ip)) 



¥=3 



(32) 
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By matching boundary conditions in the left potential we obtain two equations. The 
first one is the exact equation for the position of the levels near the barrier top in the 
case of small viscosity limit 



COS 



j d^yj2M {E - U{^)) - J dip^2M {E - U{^)) j + (1 + exp(-7rA))^/^ • 

ln0^+ j d^y/2M{E-U{ip)) + j d^^y2M{E-U{ip))] =0 



A A 
' cos I X + 2 + 2 



(33) 



V3 



V>1 



In eq. (l33!l phase shift x is defined by the expression 



1 + iX 



2^e(--A/4) 



Xix) 



(34) 



where T{x) is the Euler gamma function and the phase x can be presented in the form 



X 



X 



arctan 



A 



2k 



X 



2k + 1 



(35) 



In eq. fl35|) iplx) is the Euler function and 
1' 



-C-21n2 = -1.96351 



The second equation gives the value of the coefficient B as follows 



B = exp 



ttA 



'1 + exp (-ttA))^''^ sin I X + 



A A 



ip2 f>A 

In 0) + J dip^/2M{E -U{ip)) + j dip^/2M {E - 



sm 



dip^/2M{E -U{ip)) - / d^^2M{E -U{^)) 



<^3 



In order to obtain two close levels near to the barrier top (but not too close) the 
parameter A should be A > 1. In such a case, the normalization factor G can be 
approximated as 




dip 



If A, 



v/2M(E-f/(<^)) 



5^ 



dip 



Vi 



l2M{E -\J{ip)) 



(37) 



Since we consider the energy E close to the barrier top, it's possible to use expressions 
with explicit energy dependence (perturbation theory over the quantity {Uto-p — E) with 
separation of the singular terms): 



tP2 



dip^2M{E -U{^)) 
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\ \/ {Utop 

(v?4 - ftop) y/Ui \ , 1 
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In 



\/Utop — E j 2 



In eqs. (|551) . (15^ quantities ipi, ip^ are defined as 

ipi = {E = Utop) 

P>4 = ip4 {E = Utop) 

Energies E^, Eji of states can be found by using eqs. fl55]) . (1551) . fl5I?l) . provided 

tliat "^L, are referred to the barrier top. In this case wavefunctions "^l, can be 
described by using the quasiclassical approximation, so that we find 

sin j f + / dip^2M{E-U{ip)) 
Gl {2M{E-U{v))f" ^ ^ 



sin f + / dp>^2M{E-U{p>)) 

= i 1^ ^ ^ (41) 

Gr {2M{E-U{^))f/' ^ ^ 

where the factors Gl and Gr are defined as 

.2 ^ [ df 



G = — (42) 
2j J2M {E-U{ip)) 



V4 

Gl = ^ [ , (43) 
2j y/2M{E-U{^)) 



(44) 



The wave function \[^o of the ground state can be taken in the form 

(2Mf/™")^/® (" 1„ <i^^2M{u{v)-U{vrn)^ 

where [/{"g" = U{(p'^'!^) and v^™" coordinates corresponding to the position of 

the left and right minimum in the potential U{p>) 
fdU_\ _ ^ 

TTt,^\ TT(,^'min\ < TTmin/,„ ,„min\2 i 

U[<p) = U{ip^^2 ) + Ui^2 W - Vi,2 ) + ■■■ 

4. The energy spectrum in the vicinity of the crossing point. 

In order to describe the spectrum near to the crossing point we introduce two functions 
$1 and $2j defined as: 

$1 = y dy^^/2M {E - U{ip)) + ^ (^1 + In 0^ ^ = j dip ^2M {Utop - U{ip)) 



Resonance phenomena in Macroscopic Quantum Tunneling: the small viscosity limitl2 

1 a/ V^top - '^1 




21/4 



(45) 



$2 = j^' dip^2M{E-U{ip)) + ^ (^1 + In ^ = j^^' dip^2M{Utop~U{ip)) 

— [Utop — ^> ' 



2 




Now, we suppose that for some value of the external parameters (p^,Uo,Pl there is a 
point characterized by Eq, Xq = \ (Eq) (defined in eq j3T|) such that 

$1 (f/o, Pl, Eo, Ao) + (Ao) = ^ + vrfci (47) 

$2 (f/o, </^°, Eo, Ao) + (Ao) = ^ + 7ck2 (48) 

where fci, A;2 are integer numbers. The equations (46), (47) determine the "crossing 
point". In fact, we are interested to find the external parameter cpx such that the two 
energy levels Eji and Ef2 are close. As a consequence, next to this special point, we 
can use Taylor expansion for the functions $i and $2 and we can write: 

1 TT 

$1 + -X = - + Trfci + ai6^^ + (3i5E (49) 

1 TT 

$2 + 2^ = 2 + + 012^^^ + I325E (50) 
where 

E = Eo + 6E, ip^ = ipl + 6ip^ (51) 
and the quantities ai, 02, Pi, P2 can be found from eqs.(l2]), (l45|l . They are given in 



Appendix B. From eqs. (l33!) . (jUj), (!50ll we obtain the equation to calculate the energy 
spectrum for two close levels near to the barrier top 

PMSEy + a,a2{6ip^)^ 

+ {aiP2 + a2pi)S^JE - ^e(-"^) = (52) 
The solutions of this equation are two hyperboles 

6E = — 

2P1P2 

Note that Eql53] gives the position of two close levels as function of the external param- 
eter ipx- Two levels can be defined close if 6E « hflp (the plasma frequency, defined 



(ai/?2 + a2pi)5<fx ± V {oiip2 - a2piyi5ipxy + P1P2 exp(-7rA) (53) 
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Figure 3. Transition probability W vs. (px for different values of the pumping 
frequency v. W presents three peaks: the central one is due to the resonance tunneling 
and two lateral ones are due to the resonant pumping of the two levels close to the 
barrier top. Curves are obtained by using the following parameters for the numerical 
simulation: /3l=:1.75, L=210 pH, C=0.1 pF, and i?e//=8 Mfi. 

in the introduction) and this occurs if the parameter exp(— vrA) is small, that is if A is 
of the order of one. 



5. Numerical simulations and preliminary measurements 

In this section we suggest a set of SQUID parameters useful in order to observe exper- 
imentally the effect here studied with this new theoretical approach. The experiment 
can be realized by following the experimental scheme described in [8]. The external 
parameter (y^^. biasing the rf-SQUID is varied as a function of the time until a switching 
from one fluxoid state to another one is observed. The statistical distribution of the 
switching value (p^ can be measured by repeating the process many times (about 10'^ 
measurements). From the switching value distribution is straightforward to obtain the 
escape rate W and compare data and theory. The first significant comparison is then 
the dependence of the escape rate as a function of ip^, showing the characteristic peaks 
due to ELQ and MQC phenomena. For these simulations we consider a temperature 
T=50 mK and some reasonable rf-SQUID parameters, as: /5i=1.75, L=210 pH, C=0.1 
pF, and Reff=8 MQ. 

In the extremely underdamped limit, the number of peaks of the transition probability 
W vs. ipx, can be three or one depending on the pumping frequency value u. 

In the three peaks curves (FigJSj), the central peak is always connected with the re- 
sonant tunneling, while the two other peaks are connected with the resonant pumping. 
The central peak is due to the coherent tunneling between levels /i and f2- As a con- 
sequence, it's fixed in position and height, even if we change the frequency z/, since it is 
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correlated to the anticrossing point, that does not depend on the external irradiation. 
On the contrary, the other two peaks have a shift depending on the frequency. The left 
side peak corresponds to a transition from the ground state in the left well to the state 
/i, while the right side peak is connected to the transition from the ground state to the 
state /2. Of course these transitions occur at values of the external flux ip^ which depend 
on the pumping frequencv, as shown in FigJH Infact, if we represent the external field 

ve is possible only when the 

t in a small frequency range 



0,344 0.34S 34S 0,347 0.348 0.^3 

%. 

Figure 4. Curves Ej^ — Eq vs. (px and Ef^ — Eq vs. (px- AH the energies are referred 
to the barrier top. The horizontal hne is the external field of frequency v, it crosses 
the energy curves Ef^-E^ and Ef^-Eo in two points. The vertical line represents the 
coordinate corresponding to the anticrossing point. Comparison of Figs [3] -[5] shows 
that the first peak correspond to the resonant pumping from the ground state in the 
left well to the state with energy Ej-^ while the third peak corresponds to the resonant 
pumping from the ground state from the ground state to the state with energy Ef^. 
For all the parameters values are: Pl — 1.75, C = 0.1 pF, L = 210 pH. 

around the anticrossing point. For all other frequencies we find the one peak curve 
(FigIS]), where the only peak is due to the resonant tunneling. 

For these simulations we have considered pumping frequencies ranging from u = uj/2tt 
= 25.6 GHz to u = 26.8 GHz. These values assure that the excited pumped level in the 
left potential well is close to some level in right potential well. 

Moreover it is interesting to study the dependence of the peaks on the effective dissi- 
pation described in the RSJ model by the effective resistance -Re//, and this is shown 
in Figini As expected [27], by decreasing dissipation, that is by increasing -Re//, peaks 
resolution is enhanced. 

Also the (3l parameter can be varied in the experiments [9j and therefore we present 
some numerical simulations to study the transition probability W vs.ip^ for two different 
Pl values (see figd figlDand fig|9]). It's worth noting that to vary (3^ for the rf-SQUID, 
we vary only the value of the inductance L (see eqJl]) and a variation of Pl corresponds 
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be W presents only one peak due 
^■*-1»76 -. yy using the following parameters 
_d ::;=0.1 pF, and Reff=S MO. 
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Figure 6. Transition probability W vs. cpx for different values of the effective 
resistance Reff- The plot was obtained by using the following parameters for the 
numerical simulation: /?i=1.75, L=210 pH, C=0.1 pF, u= 25.756 GHz. 



to the variation of the height of the potential barrier. We stress that, when we increase 
the Pl value, the transition probability becomes smaller and we need to increase the 
external parameter to observe the peaks. 



Finally we also made a study of the dependence of the peaks distribution and behavior 
on the capacitance C. We observed that by decreasing the C value, the peaks in the 
transition probability become higher and more enhanced. 

The best set of values for these parameters is when the (3l and the capacitance C values 
are not too large. 

The predicted process can be observed by escape rate measurements. Since we want 
study Macroscopic Quantum Phenomena, we projected and realized devices with a low 
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Figure 7. Transition probability W vs. for /3l=2.15, L=258 pH, C=0.1 pF, 
Reff=8 Mfl and z/=32.818 GHz. 
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Figure 10. Dots are experimental data for Reff in our devices at different 
temperatures, while the line is the linear fit, corresponding to the expression Reff = 
Ro exp(A/A:Br). 



dissipation level. The preliminary measurements on our devices have confirmed a really 
low dissipation level [5], with a really high effective resistance, as shown in FiglTOl These 
results indicate that our devices are good candidates for quantum measurements. 

6. Conclusions 

The results presented here describe the small viscosity limit in the problem of resonant 
quantum tunneling for an rf-SQUID, whose parameters can be determined with 
independent measurements. All the physical quantities here considered are very sensitive 
to the external parameters. 

It's worth noting that we considered levels below the top of the classical energy barrier. 
This condition is essential, since in this case the levels (in the absence of coherence) can 
be exclusively associated with one well and moreover it's a necessary condition in order 
to observe tunneling phenomena. 

The transition probability function W from the left potential well to the right one has 
been calculated by varying cp^ and for different values of the pumping frequency u, of 
the parameter and for different values of the effective resistance Reff- The results 
obtained here suggest that the resonant quantum tunneling is a convenient tool for 
investigating macroscopic quantum coherence phenomena. 

The dependence of the transition probability W on the external parameter shows 
three peaks (Figj3]). The first is connected with the resonant tunneling and the others 
two are associated with the resonant pumping toward the two non-localized energy levels. 
The relative positions of these two peaks strictly depend on the pumping frequency and 
on the external parameter biasing the rf-SQUID. Moreover the transition probability W 
is obtained as a function of the potential parameters and of the shunt resistance Reff- 
Preliminary experimental measurements on our devices confirmed a low dissipation 
level [5], with a high effective resistance, as shown in FigJTOl These results, albeit 
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still unconfirmed by a wider and more accurate experimental campaign, would suggest 
potential for interesting observations regarding the resonant phenomena in the presence 
of an external microwave irradiation of proper frequency. 
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Appendix A. 

From eqs.(I7j), (|T3l) we obtain two equations depending on the quantities Vi and J^i 



^1 



Ef2 - E 



h 



h 



Ef2 - Ef^ 

h 



16^2 
16e2 



{0Wi){f2\mwi 





/2 



- 7172 + Wl)Wl\ + z(7i + 72) - 
(0|v^|/i)(/2|^|0)W^;,°^, 



- 7172 + Wl^f.Wl^ - z(7i + 72) [uJ 



Solution of equations (lA.ip is 

^ (0|^|/l)(/2|^|0) 

^1 



- 47172 - 22(71 + 72) 



Ef2-Efi 



ji j2 h h 



Ef2-Efi 



-2z72 



Ji 72 



"-^)" - 7172 + W^^lWll + z(7i + 72) - ^ 



(A.3) 
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UJ — 



Eh-Eo\ 



ZJ2 



j - 7172 + /,W^r/i - ^(71 + 72) [UJ 
(0|v^|/i)(/2|^|0) 



16e2 fEj^-Ef^ 



- 47172 - 2^(71 + 72)^^^^ + WllllWllll 



Ef,-Eo 



-2^71 



^^f2 h 



- - ^) - 7172 + Wl\Wl\ - Kli + 72) - ^ 

</>r/2 



(A.4) 



UJ 



Ef2-Eoy 



-7i72 + W^r/2^/2/i+^(7i + 72)(c^-^ 

The exphcit expression for the transition matrix elements ^j^j^ s^nd W^/^/j is given 
below 



w/1/2 _ 



IT 



w 







1 + tanh 
1 + tanh 



Ef2 - Ef, 
2kBT 

Ef2 - Ef^ 



Ef^ — E 



-coth 



~ coth 
27r 



2kBT 
Ef2 - E 



fi 



(0|exp^2 |0)(/i|exp-2 \f,) + (0|exp-*2 |0)(/i|exp^2 j/^) 
1% = -P 



^/2 - ^/l 



2A;ijT 



'^'^/l/2 



(A.6) 
(A.7) 
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Appendix B. 

For the potential, given by eq.(l2]) we have 
V^top = <^x + Pl sin(v9top) 

dipx 
dU 



I3l cos{iptop) - 1 



f^l = Y COs{(ptop) - 1) 

d(f^ 2 Pl cos(v?top) - 1 



Utop — Uq 



da 



top 



-^{^top - ^xf + h COS iftop 

(3l sin LP top - {ptop - Px) 



- {Ptop - <Px 



(3l cos(v?top) - 1 

From equations fl29l) . fl32l) we obtain all other derivatives that should be found for 
calculation of the quantities ai,2, Pi,2- 

'dE 
dX 



2M 



d(px 



Ur 



2M 



_Pl COs{(ptop) - 1 
oo 



13 L sin ^top ( 1 - ^^^^ 



E 



{Ptop - Vx) 



k=0 



{2k + l)i{2k + ly + X^) 



(B.l) 

{iptop - ^•^) 

(B.3) 



Now by using equations (H5l) . (H6l) . (IB. II) . ( IB. Ill we obtain the value of all coefficients 

"l, "2, A, P2- 

ft op 



Pi — ^TTT + 




M 



+ 



dE 2 dX dE 

{MUiy/\ptop-Pi) 



dip 



\l 'Ptop - Vi 



\/Utop - U {p) {p>top - p) \/Ui{p - pi) 



(B.4) 
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Appendix C. 

Transition matrix elements between states close to the barrier top can be calculated in 
quasiclassical approximation if the parameter A (see eq. fl28l) ) is A > 1. Consider first 
the matrix element (L|C|/i). With the help of eqs.(l28]), (HQD, (HI]) we obtain 



^2 



(Ef,) sin 7r/4+ / d^^{2M{EL-U{ip))) 



( 



(^ICI/i) 



GlGe 



fi 



d(f- 



(2M {El - U{ip))f" (2M {E^^ - U{ip))) 



1/4 



Csin 



71/4:+ j d^^J{2M{Ef 



that is equal to 



(^ICI/i) 



dip 



■.Ccosl{Ef,-EL) r 



fi{Ef^) 

Mdip 



. ^{2M{Ef,L-U{v))) V 'JMEf,.)V^M{EjJ^UM, 

To calculate the last integral we can use the time variable t according to the classical 
equation of motion 

As a consequence, eq. (IC.2p becomes 

Ti/2 

(^ICI/i) = TTTTT^^ f dtCcos(^^te] (C.4) 



^ 

where Ti is the period of the classical motion in the left potential well with energy E 
and i is the number of states between energy levels Ef-^ and E^ plus one. In our case 
i = 1. The initial value for (p is f{0) = fi- 

dcp 



Ti = 2M 



V(2M(E-f/(^))) 



(C.5) 



In the same way we obtain transition the matrix element {R\(\fi). 



In eq. OC.ip ipi{Ef^) is the first crossing point of the energy level Ef^ in the left well 
and (^2(-E'/i) is the second crossing point of the energy level Ef^ in the left well. 
In the same way in eq. flC.2p ipi^Ef^i) is the first crossing point of the energy level Ef^L 
in the left well, where Ef^i is defined as 

-E'/i + El 



2 

Consider now the transition matrix element (/2I Cl/i)- From eq. fl28|) we find 

r 

dip- 



(C.6) 



(/2ICI/1) 



2Gf^Gf2 



ip4 

+ Bf,Bf, j dip 



c 



^{2M{E-U{ip))) 



(C.7) 



(C.2) 
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( E -\-E 1 

where (/'i, 2,3,4 are the "turning points" and E here is E = ^ ^ 
By using eq. (]C.3|) we reduce the expression (]C.7p to the form 

'Ti/2 T2/2 

j dtC + Bf-^Bf^ J dtc 



(/2I Cl/i 



2MGf,Gf, 



(C.8) 



In eq. flC.81) the first integral is taken over the left potential well and the second one 
over the right potential well and Ti, T2 are the periods of the classical motion in the left 
and in the right well respectively. 

To improve eq. flC.8l) we should take into account ortogonality of wavefunctions relative 
to the energies Ef^^, Ef^ . As result we obtain 



(/2I Cl.f1) = ft 



IMG f-^ 



T2 
T1+T2 



Ti/2 



dtC- 



T1+T2 



T2/2 



dtc 



(C.9) 



Acknowledgments 

The research of one of us (Yu.N.O.) is supported by Cariplo Foundation-Italy, CRDF 
USA under Grant No. RP1-2565-MO-03 and the Russian Foundation of Basic Research. 
This work has been partially supported by MIUR under Project "Reti di Giunzioni 
Josephson per la Computazione e I'lnformazione Quantistica-JOSNET" . 



References 

[1] A. O. Caldeira and A. J. Leggett. Phys. Rev. Lett, 46:211, 1981. 

[2] D. B. Schwartz, B. Sen, C. N. Archie, and J. E. Lukens. Phys. Rev. Lett., 55:1547, 1985. 
[3] S. Washburn, R. A. Webb, R. F. Voss, and S. M. Paris. Phys. Rev. Lett., 54:2712, 1985. 
[4] R. Rouse, S. Han, and J.E. Lukens. Phys. Rev. Lett., 75:1614, 1995. 

[5] V. Corato, S.Rombetto, C. Granata, E. Esposito, L. Longobardi, M. Russo, R. Russo, B. Ruggiero, 

and P. Silvestrini. Phys. Rev., B70:172502, 2004. 
[6] J. M. Martinis, M. H. Devoret, and J. Clarke. Phys. Rev. Lett., 55:1543, 1985. 
[7] J. M. Martinis, M. H. Devoret, and J. Clarke. Phys. Rev., 35:4682, 1987. 
[8] P. Silvestrini, V. G. Palmieri, B. Ruggiero, and M. Russo. Phys. Rev. Lett, 79:3046, 1997. 
[9] J. R. Priedman, V. Patel, W. Chen, S. K. Tolpygo, and J. E.Lukens. Nature (London), 406:43, 

2000. 

[10] J. B. Majer, P. G. Paauw, A.C. J. ter Haar, C. J. P.M. Harmans, and J. E. Mooij. Phys. Rev. 
Lett, 94:090501, 2005. 

[11] J. Claudon, P. Balestro, P. W. J. Hekking, and O. Buisson. Phys. Rev. Lett, 93:187003, 2004. 
[12] M. H. Devoret and J. M. Martinis. Quantum Entanglement and Information Processing. - Les 

Houches Summer School Proceedings 79. Elsevier, San Diego, CA,, 2004. 
[13] N. Groenbech-Jensen and M. Cirillo. Phys. Rev., B70:214507, 2004. 
[14] A.I.Larkin and Yu.N.Ovchinnikov. J. of Low Temp. Phys., 63:317, 1986. 
[15] Yu.N.Ovchinnikov, P. Silvestrini, V.Corato, and S.Rombetto. Phys. Rev., B71:024529, 2005. 
[16] Yu.N.Ovchinnikov and A.Schmid. Phys. Rev., B50:6332, 1994. 
[17] D. Averin. Solid State Commun., 105:659, 1998. 
[18] D. Averin. Phys. of low temp. Phys., 118:781, 2000. 

[19] D. V. Averin, J. R. Friedman, and J. E. Lukens. Phys. Rev., B62:11802, 2000. 
[20] A.I.Larkin and Yu.N.Ovchinnikov. Sov. Phys. JETP, 64:185, 1986. 



Resonance phenomena in Macroscopic Quantum Tunneling: the small viscosity limit2A 



[21] A.I.Larkin and Yu.N.Ovchinnikov. Phys. Rev., B50:6332, 1994. 
[22] U. Fano. Rev. of Mod. Phys., 29:74, 1957. 

[23] A. Barone and G. Paterno. Physics and Applications of the Josephson Effect. Wiley, New York, 
1982. 

[24] A.I.Larkin and Yu.N.Ovchinnikov. Sov. Phys. JETP, 60:1060, 1984. 

[25] 1. M. Ryzhik and I. S. Gradstein. Table of Integrals, Series and Products. Academic Press, New 
York, 1965. 

[26] L. Landau and L. Lifshitz. Quantum Mechanics: non-relativistic theory. Butterworth-Heinemann, 
1981. 

[27] P. Silvestrini, B. Ruggiero, and Yu.N. Ovchinnikov. Phys. Rev., B54:1246, 1996. 



